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One of the main scientific objectives of the recent Fermi mission is investigating the origin of the unidentified 
gamma-ray sources (UGSs). Despite the large improvements of Fermi in the gamma-ray source localization 
with respect to the past gamma-ray missions, about 1/3 of the gamma-ray objects detected do not have yet 
an assigned counterpart a low energies. We recently developed a new association method to identify if there 
is a 7-ray blazar candidate within the positional uncertainty region of a generic 7-ray source. This method is 
entirely based on the discovery that blazars, the largest known class of gamma-ray sources, can be recognized 
and separated from other extragalactic sources on the basis of their infrared colors. Here we summarize the 
results obtained by applying our association procedure to the unidentified 7-ray sources (UGSs) and to the 
active galaxies of uncertain type (AGUs) listed in the Second Fermi Large Area Telescope (LAT) catalog. 



1. Intorduction 

Unveiling the natu re of the Unident ified Gamma- 
ray Sources (UGSs) jAbdo et al. 2009j| is one of the 
biggest challenges in contemporary gamma-ray as- 
tronomy. In particular, according to the Second 
Fermi Large Are a Telescope (LAT) catalog [2FGL; 
iNolan et al. 20121 ] , ~ 1/3 of the 7-ray detected sources 
are still unassociated with their low energy counter- 
parts. A large fraction of the UGSs are likely to 
be blazars, the rarest class of active galaxies, be- 
cause their emission dominates the 7-r ay sky [e.g., 
iMukheriee et al. 1997tlAbdo et al. 2010j . 

However, due to the incompleteness of the cur- 
rent radio and X-ray surveys on the basis of the 
usual 7-ray association method is not always possi- 
ble to find the blazar- like counterpart of an UGS. 
Additional attempts have also been recently devel- 
oped to associate or to characterize t he UGSs using 
either pointed Swift obse rvations [e.g.. iMirabal 2 009. 
IMirabal fc Halpern 20091 or statistical approa ches 



Blazars come in two flavors: the BL Lac ob- 
jects, characterized by featureless optical spectra 
and lower luminosity with respect to the second 
class composed of flat-spectrum r adio quasars show- 
ing quasar-like o ptical spectra jStickel et al. 1991. 
Stoke et al. 199l| . In the following we label the 
BL Lac objects as BZBs and the flat-spectrum ra- 
dio quasars as BZQs, following the nomenclature of 
the Multi frequency Catalogue of Blazars [ROMA- 
BZCAT. IMassaro et al. 20091 iMassaro et al. 20101 
iMassaro et al. 2011a [. 

On the basis of the preliminary data release 
of the Wide-field Infr ared Survey Explorer [WISE, 
see IWright et al. 2010L for more details!, we discov- 
ered that in the 3-dimensional IR color space 7-ray 
emitting blazars lie in a distinct region, well sep- 
arated from other extragalactic sources whose IR 
emission is dominated by thermal radiatio n [e.g. 



Ma ssaro et al. 201 lbl iD'Abrusco et al. 20121 . Ac 



[e.g. IMirabal et al. 2010l lAckermann et al. 2012] . 

Blazar emission is characterized by high and vari- 
able polarization, apparent superluminal motions, and 
high luminosities, coupled with a flat radio spectrum 
that steepens toward the infrared-optical bands and 
toget her with rapid flux var iability at all frequencies 
[e.g.. lUrry fc Padovani 19951 ]. Their broad band spec- 
tral energy distributions show two main components: 
the low-energy one with power peaking in the range 
from the IR to the X-ray band, and the high-energy 
show ing its maximum in t he MeV - TeV energy range 
[e.g.. lGiommi et al. 2005j . 



cording to D'Abrusco et al. (2013) we refer to the 3- 
dimensional region occupied by 7-ray emitting blazars 
as the locus, to its 2-dimensional projection in the 
[3.4]-[4.6]-[12] /im color-color diagram still maintain 
its historical definition as the WISE Gamma-ray Strip. 

This WISE analysis led to the development 
of a new association method to recognize 7- 
ray blazar candidates for t he unidentified 7-ra; 



sources listed in th e 2FGL [M assa ro et al. 2012 
IMassaro et al. 2012bl. as well as in the 4 th INTE- 
GRAL catalog [Massaro et al. 2012cj . 

Here we present the results achieved by applying 
a more conservative approach and several improve- 
ments recently made on the association procedure [see 



eConf C121028 



4 th Fermi Symposium : Monterey, CA : 28 Oct- 2 Nov 2012 



iD'Abrusco et al. 20131 for more details] , mostly based 
on the availability of the WI SE full archive 1 , avail- 
able since March 2012 [see also lCutri et al. 20121 ]. We 
analyzed all the UGSs listed in the 2FGL as well 
as the sample of the active galactic nuclei of uncer- 
tain type (AG Us) that have still unclear classification 
Nolan et al. 2012^ . We also performed an extensive 



literature search looking for multifrequency informa- 
tion on the 7-ray blazar candidates selected on the 
basis of their WISE colors to confirm their nature. 

We remark that a detailed description of our mul- 
tifrequency analysis and of the association procedure 
adopted, together with the lists of j-ray blazar candi- 
dates, is presented in details in Massaro et al. (2013) 
and D'Abrusco et al. (2013), while here we only high- 
light the major results achieved. 



2. RESULTS 

A new association method has been recently de- 
veloped based on the discovery that 7-ray emit- 
ting blazars lie in a distinct region in the WISE 
3-dimensional color space, separated from that oc- 
cupied by other extr a galactic and galactic sources 



Massaro et al. 2011bl iD'Abrusco et al. 2012j . Ac 



cording to D'Abrusco et al. (2013) the 3-dimensional 
region occupied by 7-ray emitting blazars is the locus; 
its 2-dimensional projection in the [3.4]-[4.6]-[12] /mi 
parameter space, retains its historical definit ion of 
WISE Gamma-ray Strip [Massaro et al. 2011blj . Ad- 
ditional improvements, mostly based on the WISE 
all-sky data rele ase, available since March 2012 
Cutri et al. 2012j . and on a new parametrization 



of the locus in the parameter space of its princi- 
pal components have been subsequently developed 
[D'Abrusco et al. 2013| . 

Then we applied our new association procedure 
searching for new 7-ray blazar candidates in the two 
samples: the unidentified gamma-ray sources (UGSs), 
and the active gala xies of uncertain ty pe (AGUs), as 
listed in the 2FGL [Nolan et al. 2012j ]. We also per- 
form an extensive archival search to see if the sources 
associated with our method, show additional blazar- 
like characteristics; as, for example, the presence of a 
radio counterpart and/or of a spectrum that could be 
featureless as for BZBs or similar to those of broad- 
line quasars as generally occurs in BZQs. 

We found 62 7-ray blazar candidates for the UGSs 
without any 7-ray analysis flag and 49 for those with 
7-ray analysis flag, out of a total of 590 sources in- 
vestigated. For the AGU sample, we confirmed the 
blazar-like nature of 87 out 210 of AGUs analyzed 
on the basis of their IR colors. Figure [T] shows the 



3-dimensional color plot comparing the IR colors of 
the selected 7-ray blazar candidates with the blazar 
population that constitutes the locus. 
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Figure 1: The 3D representation of the locus (known 
7-ray blazars are indicated in yellow) in comparison with 
the selected 7-ray blazar candidates: UGSs (red) and 
AGUs (black). 

We searched in the following major radio, infrared, 
optical and X-ray surveys as well as in the NASA 
Extragalactic Database (NED) 2 for possible counter- 
parts of our 7-ray blazar candidates, selected with the 
WISE association method, to verify if additional in- 
formation could confirm their blazar-like nature. For 
the radio count erparts we used the NRAO VLA Sky 
Survey [NVSS; ICondon et al. 19981 - N], the VLA 
Faint Im ages of the Radio Sky at Twenty- Centim eters 
[FIRST; iBecker et al. 199,1 White et al. 19971 - F], 
the Sydney Univer sity Molonglo Sky Survey [SUMSS; 
iMauch et al. 20031 - S] and the Australia T elescope 
20 GHz Survey [AT20G: lMu"rphv et al. 20H - A]; for 
the infra red we used the Two Micron All Sky Survey 
[2MASS: ISkr"utskie et al. 200rl - M] since each WISE 
source is already associated with t he closest 2MASS 
source by the default catalog [see ICutri et al. 20121 
for more details]. Then, we also searched for 
optical counterparts, with possible spectra avail- 
able, in the Sloan Digital Sky Survey [SPSS; e .g. 
lAdelman-McCarthv et al. 20081 iParis et al. 20121 - 
s], in the Six-degree-Field Galaxy Redshift Sur- 



1 http: / /wise2. ipac.caltech.edu/docs/release/allsky/ 



2 http:/ /ned. ipac.caltech.edu/ 
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vey [6dFGS; iJones et al. 2004 iJones et al. 20091 - 
6]; while for the high energy we looked in the 
soft X-rays using the ROSAT all-sky survey [RASS; 
IVoges et al. 19991 - X]. We also s earched in the 
USNO-B Catalog (Monet et al. 2003j for the optical 
counterparts of our 7-ray blazar candidates 

A significant fraction (i.e., ~ 36%) of the WISE 
sources associated with our method with UGSs have 
a radio counterpart, more than 50% are also detected 
in the 2MASS catalog as generally occurs for blazars, 
and more than ~10% appear to be variable accord - 
ing to the WISE analysis flags Cutri et al. 2012} . 
Notably, more than 90% sources for which an opti- 
cal spectrum was available in literature clearly show 
blazar-like features, being either featureless or having 
broad emission lines typical of quasars. As generally 
expected for 7-ray blazars a handful of the selected 
candidates are also detected in the X-rays. A deeper 
investigation of their X-ray cou nterparts will be ad- 
dressed in a forthcoming paper [Paggi et al. 2013j . 

Our results are in good agreement with those based 
on different statistical approaches like the Classi- 
fication Tree and the Logistic regression analyses 
[Ackermann et al. 2012j |. In particular, 23 out of 28 
UGSs that we associate to a 7-ray blazar candidate are 
also classified as active galaxies by the above methods 
at high level of confidence. 

Finally, we emphasize that additional investigations 
of different samples of active galactic nuclei, such as 
Seyfert galaxies, are necessary to study the problem 
of the contamination of our association method by 
extragalactic sources with infrared colors similar to 
those of 7-ray blazars. Moreover extensive ground- 
based spectroscopic follow up observations in the op- 
tical and in the near IR would be ideal to verify the 
nature of the selected WISE sources and to estimate 
the fraction of non-blazar objects, similar to the recent 
studies perf ormed for the unidentified INTEGRAL 
sources [e.g.. IMasetti et al. 20081 IMasetti et al. 20091 



IMasetti et al. 20101 IMasetti et al. 2012} . 
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